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The microscopic description of the spin-Peierls transition in pure and doped CuGeOs is developed 
taking into account realistic details of crystal structure. It it shown that the presence of side-groups 
(here Ge) strongly influences superexchange along Cu-O-Cu path, making it antiferromagnetic. 
Nearest-neighbour and next-nearest neighbour exchange constants J„„ and Jnnn are calculated. 
Si doping effectively segments the Cu02-chains leading to J„„(Si) ~ or even slightly ferromag- 
netic. Strong sensitivity of the exchange constants to Cu-O-Cu and (Cu-O-Cu)-Ge angles may 
be responsible for the spin-Peierls transition itself ("bond-bending mechanism" of the transition). 
The nature of excitations in the isolated and coupled spin-Peierls chains is studied and it is shown 
that topological excitations (solitons) play crucial role. Such solitons appear in particular in doped 
systems (Cui-^ZniGeOs, CuGei-iSi^Os) which can explain the Tsp{x) phase diagram. 
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I. INTRODUCTION 

The low-dimensional materials are known to be very 
susceptible to various instabilities, such as formation of 
charge- or spin-density waves. Probably, the first one dis- 
cussed is the famous Peierls instability of one-dimensional 
metals: lattice distortion with the new lattice period 
2tt/Q, where the wave vector Q — 2kp (if there is one 
electron per site, the lattice dimerizes). The lattice dis- 
tortion opens a gap in the electron spectrum at the Fermi 
surface, so that the energies of all occupied electron 
states decrease, which drives the transition. It is also 
known that this instability survives when we include the 
strong on-site Coulomb repulsion between electrons (the, 
so called, Peierls- Hubbard model). 
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Here the first term describes the dependence of the elec- 
tron hopping integral U^i+i on the change of the distance 
ui—ui+i between the neighbouring ions and the last term 
is the lattice energy (which after quantization becomes 
'^q^qb^bq). Thc dimcnsionlcss electron-lattice coupling 
constant A — Aa'^ / {irtoK) determines the magnitude of 
the lattice distortion and the energy gap. 

When Coulomb repulsion is strong, U ^ t^, and there 
is one electron per site, we are in the limit of localized 
electrons (Mott-Hubbard insulator) with effective anti- 
ferromagnetic interaction (Spin-Peierls model). 
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where the exchange constant J/.;+i — Jq + a'{ui — u;+i), 
Jo = 4tQ/C/ and a' — 8toa/U. The dependence of J/.;+i 
on the distance between neighbouring spins again leads to 
an instability, as the result of which the spin chain dimer- 
izes. Physically it corresponds to a formation of singlet 
dimers — the simplest configuration in the valence bond 
picture. This transition, known as the spin-Peierls (SP) 
transition, was extensively studied theoretically and 
was previously observed experimentally in a number of 
quasi-one-dimensional organic compounds, such as TTF- 
CuBDT {Tsp = 12K) or TTF- AuBDT {Tsp = 2.1K) 

i- 

Recently the first inorganic spin-Peierls material 
CuGeOa was discovered . Since then much experimen- 
tal data on this material, both pure and doped (mostly 
by Zn and Si), was obtained. The spin chains in this 
compound are formed by CUO4 plaquettes with common 
edge (See Fig.l). They apparently play the main role in 
the spin-Peierls transition with Tsp = UK. However, as 
we will discuss below, the interchain interaction is also 
very important here. The interchain coupling is provided 
both by Ge ions (along 6-axis of thc crystal) and by the 
well separated apex oxygens (along a-axis; direction of 
the chains coincides with the c-axis of a crystal). 

Experimentally it is established that the strongest 
anomalies in CuGeOa, both in the normal phase and at 
the spin-Peierls transition, occur not along the c-axis, 
but, rather unexpectedly, along the other two directions, 
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the strongest one found along the 6-axis For in- 

stance, the anomahes in the thermal expansion coelRcient 
and in the magnetostriction along the 6-axis are several 
times stronger than along the direction of the chains . 

Further interesting information is obtained in the 
studies of doping dependence of various properties of 
CuGeOa. It was shown that the substitution of Cu by 
nonmagnetic Zn, as well as Ge by Si, leads initially to a 
rather strong reduction oiTsp 10 1, which according to 
some recent data Q is flattened out at higher doping 
level. Simultaneously, antiferromagnetic order develops 
at lower temperatures, often coexisting with the SP dis- 
tortion 

The aim of the present investigation is to provide mi- 
croscopic picture of the properties of CuGeOa taking into 
account realistic details of its structure. We will address 
below several important issues: 

• The detailed description of the exchange interac- 
tion (Why is nearest neighbour exchange antiferro- 
magnetic?); 

• The sensitivity of the exchange constants to dif- 
ferent types of distortion and the resulting from 
that microscopic picture of the spin-Peierls tran- 
sition, which may be called bond-bending model 
(Why the anomalies are strongest in the perpen- 
dicular directions? Why is SP transition observed 
in CuGeOa and not in many other known quasi- 
one-dimensional magnets?); 

• The nature of elementary excitations in spin-Peierls 
systems in general (Are they the ordinary singlet- 
triplet excitations? How are they influenced by the 
interchain interaction?); 

• The mechanism by which doping affects the prop- 
erties of SP system (Why is the effect so strong? 
Why does the system develops antiferromagnetic 
order upon doping?). 

These questions are raised by the experimental obser- 
vations, and we hope that their clarification will help 
both to elucidate some general features of SP transitions 
and to build the detailed picture of this transition in 
CuGeOa. 



II. EXCHANGE INTERACTION IN CUGEO3. 
ROLE OF SIDE-GROUPS IN SUPEREXCHANGE 

The first question we would like to address is: why 
is the nearest-neighbour Cu-Cu exchange interaction an- 
tiferromagnetic at all? The well-known Goodenough- 
Kanamori-Anderson rules state, in particular, that the 
90°-exchange between two half-filled orbitals is ferromag- 
netic. In CuGeOs the Cu-O-Cu angle 6 in the superex- 
change path is 98°, which is rather close to 90°. Usu- 
ally, the antiferromagnetic character of the exchange is 



attributed to this small 8° difference. Our calculation 
[13], however, shows that it is not enough: even in this 
realistic geometry the exchange constant for the near- 
est neighbour (nn) spins for realistic values of parame- 
ters (such as copper-oxygen overlap, magnitude of the 
Coulomb interaction on copper and oxygen, Hund's rule 
intraatomic exchange, etc) is still slightly ferromagnetic, 
J^,^ = -0.6meV. 

To explain the observed values of the idea was put 
forward in p^ ] that the antiferromagnetic coupling may 
be enhanced by the, initially ignored, side-groups effect 
(here Ge). As it is clear from Fig.l, there is a Ge ion 
attached to each oxygen in Cu02 chain. The Coulomb 
interaction with Ge^+ and the hybridization of 2py or- 
bital of oxygen with Ge (see Fig. 2) destroys the equiva- 
lence of p2;-orbitals (shaded) and p^-orbitals (empty) of 
oxygen, which for 90°-exchange was responsible for the 
cancellation of the corresponding antiferromagnetic con- 
tributions to superexchange. As a result, the exchange 
with partial delocalization into Ge may become antiferro- 
magnetic even for 90° geometry (for a detailed discussion 
see ) . The calculation gives a reasonable value for the 
nearest- neighbour exchange interaction: = 11.6meV 
(the experimental value is 9 15meV, depending on the 
procedure of extraction of J,'^„ from the experimental 
data) . 

We also calculated other exchange constants using a 
similar approach. For the interchain interaction along 
b and a axes we obtained J,^„ — 0.7meV and J°„ = 
-3 • lO^^nieV, so that J^/J^i^ « 0.06, and J^JJ^i„ ~ 
—3 • 10 ^. The experimental values are: JnnlJnn ~ 0-1 
, JnnlJnn ~ —0.01. Thus our theoretical results are not 
so far from the experiment for the interchain interaction 
in the 6-direction and too small for a-axis. We note, how- 
ever, that the ferromagnetic exchange in a-direction is in 
any event very weak and a small variation of parameters 
used in our calculation can easily changes this value quite 
significantly. 

More interesting is the situation with the next-nearest- 
neighbour (nnn) interaction in the chain direction J^„„. 
As is clear from Figs.l and 2, there is a relatively large 
overlap of the p,^. orbitals on neighbouring plaquettes, 
which leads to a rather strong antiferromagnetic nnn cou- 
pling. Our calculation gives 7 = J^„„ /J^ « 0.23 0.3. 
From the fit of x{T) curve Castilla et al |1J] obtained the 
value 7 w 0.25. Note also that a sufhciently strong nnn 
interaction may lead to a singlet formation and creation 
of a spin gap even without the spin-lattice interaction. 
Such a state is an exact ground state at the Majumdar- 
Ghosh point 7 — 0.5 [16]. The critical value for appear- 
ance of a spin gap is 7 w 0.25 [Q. Thus, from both the 
fit to experimental data and our calculations it appears 
that CuGeOs is rather close to the critical point, so that 
one can conclude that both the frustrating nnn interac- 
tion and the spin-lattice interaction combine to explain 
the observed properties of CuGeOs (see also [H). 
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Anticipating the discussion below, we consider here the 
modification of the exchange constants caused by doping. 
In particular, we calculated the change of J^„ when Ge 
ion attached to a bridging oxygen is substituted by Si. As 
Si is smaller than Ge, one can expect two consequences. 
First, it will pull closer the nearby chain oxygen, some- 
what reducing the corresponding Cu-O-Cu angle 9. The 
second effect is the reduced hybridization of 2py orbital 
of this oxygen with Si. According to the above consid- 
erations (see also [^) both these factors would dimin- 
ish the antiferromagnetic nn exchange. Our calculation 
shows that for realistic values of parameters the re- 
sulting exchange interaction becomes either very small 
or even weakly ferromagnetic, = ± ImeV. Thus Si 
doping effectively interrupts the chains similar the effect 
of substituting Cu by Zn. This result will be used later 
in section 5. 

III. BOND-BENDING MODEL OF THE 
SPIN-PEIERLS TRANSITION IN CUGEO3 

We return to the discussion of the exchange interaction 
and its dependence on the details of crystal structure 
of CuGeOa. As follows from the previous section, the 
largest exchange constant J^^^ is very sensitive to both 
Cu-O-Cu angle 6 and to the side group (here Ge). As 
to the second factor, one has to take into account that, 
contrary to a simple model of CuGeOa shown in Fig. 2, 
in the real crystal structure Ge ion lies not exactly in 
the plane of Cu02 chain: the angle a between Ge and 
this plane is ^ 160°. The actual crystal structure may 
be schematically depicted as in Fig. 3, where the dashed 
lines represent Cu02-chains. 

One can easily understand that is also very sensi- 
tive to a Ge-Cu02 angle a. The influence of Ge, which 
according to the above consideration gives an antiferro- 
magnetic tendency, is the largest when a = 180°: just in 
this case the inequivalence of 2px and 2py orbitals shown 
in Fig. 2, which is crucial for this effect, becomes the 
strongest. On the other hand, if, for instance, a = 90° 
(i.e. if Ge would sit exactly above the oxygen) its inter- 
action with 2px and 2py orbitals would be the same and 
the whole effect of Ge on J^„ would disappear. Thus 
bending GeO-bonds with respect to Cu02-plane would 
change (it becomes smaller when a decreases). 

These simple considerations immediately allows one to 
understand many, at first glance, strange properties of 
CuGeOa mentioned in the introduction JlSf . Thus, e.g. 
the compression of CuGeOa along the 5-direction would 
occur predominantly by way of decreasing of Ge-(Cu02) 
angle a, while the tethrahedral 0-Ge-O angle (j) is known 
to be quite rigid. Such a "hinge" or "scharnier" model ex- 
plains why the main lattice anomalies are observed along 
the 6-axis and why the longitudinal mode parallel to b 
is especially soft . Within this model one can also natu- 



rally explain (even quantitatively) the fact that the mag- 
netostriction is also strongest in the b-direction j|] . If we 
assume that the main changes in the lattice parameters 
occur only due to bond bending {i.e. due to the change 
of angles, while bond lengths remain fixed), we obtain 
the following result for the uniaxial pressure dependence 
of J = J^„ (l|: SJ/SPb = -1.5mcV/GPa, which is close 
to the experimental result 6J/SPb — — 1.7meV/GPa 
We can also explain reasonably well the change of the 
exchange coupling for other directions. 

This picture can be also used to explain the spin-Peierls 
transition itself. What occurs below Tsp, is mostly the 
change of bond angles ("bond-bending"), which alter- 
nates along the chains. Experimentally it was found |l|] 
that the dimerization is accompanied by the alternation 
of Cu-O-Cu angles 9. In our model J is also sensitive to 
Ge-Cu02) angle a and we speculated in Ref.l4 that this 
angle, most probably, also alternates in the spin-Peierls 
phase. Recently this alternation was observed |l7|| . 

Consequently we have a rather coherent picture of the 
microscopic changes in CuGeOa, both above and below 
Tsp'. in the first approximation we may describe the main 
lattice changes as occurring mostly due to the change of 
the "soft" bond angles. The strongest effects for T > Tsp 
are then expected along the 6-axis, which is consistent 
with the experiment. The same bond-bending distortions 
seem also to be responsible for the spin-Peierls transition 
itself, the difference with the normal phase being the al- 
ternation of the corresponding angles in the c-direction. 

The bond-bending model allows one to explain an- 
other puzzle related to spin-Peierls transitions (discussed 
already in j2j): up to now such transitions have been 
observed only in a few of the many known quasi-one- 
dimensional antiferromagnets. There might be several 
reasons for that. The first one is that the spin-Peierls 
phase in CuGeOs is, at least partially, stabilized by the 
frustrating next-nearest neighbour interaction J^„„. The 
other factor is that the spin-Peierls instability is greatly 
enhanced when the corresponding phonon mode is soft 
enough One can see it e.g. from the expression for 
Tsp i, 

Tsp = O.SX'J . 

The spin-phonon coupling constant is 

,2 ,2 
, a a 

^ Ik ^ jmojI ' 



where = yK/M is the typical phonon frequency. 

There is, usually, a competition between the 3(i mag- 
netic ordering and the spin-Peierls phase. Apparently, 
in most quasi-one-dimensional compounds the 3d mag- 
netic ordering wins, and for the spin-Peierls transition 
to be realized a strong spin-lattice coupling, i.e. rather 
soft phonons, is necessary. Such soft phonon modes are 
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known to exist in the organic spin-Peierls compounds . 
In CuGeOa it can be rather soft bond-bending phonons, 
especially the ones parallel to the 6-axis, which help to 
stabilize the spin-Peierls phase relative to the 3d anti- 
ferromagnetic one. Nevertheless, a relatively small dop- 
ing is sufficient to make the antiferromagnetic state more 
favourable, although some other factors are also very im- 
portant here, as will become clear in the next section. 



IV. SOLITONS AND STRINGS IN SPIN-PEIERLS 
SYSTEMS 

Let us turn now to the second group of problems re- 
lated to SP systems, namely, the nature of elementary 
excitations. In the simple picture mentioned in the in- 
troduction (and valid in the strong coupling limit) the SP 
state consists of isolated dimers. For the rigid dimerized 
lattice an excited state is a triplet localized on one of the 
dimers and separated from the ground state by an energy 
gap J. The interaction between the neighbouring dimers 
gives a certain dispersion to this excitation, transforming 
it into an object similar to a usual magnon. 

If, however, the lattice is allowed to adjust to a spin 
flip, the localized triplet decays into a pair of topo- 
logical excitations. Such excitations (solitons or kinks) 
are known to be the lowest energy excitations in elec- 
tronic Peierls insulators The same is also true 
for spin-Peierls systems. Indeed, there exist two degen- 
erate ground states in a SP chain: one with singlets 
formed on sites . . . (12) (34) (56). . . , and another of the 
type . . . (23) (45) (67). . . . One can characterize them by 
the phase of the order parameter so that (f)n = in 
the first state and 0^ = tt in the second. The soliton is 
an excited state, in which the order parameter interpo- 
lates from to TT or vice versa. In the strong coupling 
limit such a state looks like . . . (12)(34)t(67)(89). . . , see 
Fig. 4. Actually, the soliton has a finite width, which (as 
the correlation length in the BCS theory) has a form. 
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Here the Fermi velocity vp ^ Ja/h is the velocity of the 
spinless Jordan- Wigner fermions, in terms of which the 
Hamiltonian (2) has a form similar to the Hamiltonian of 
the electronic Peierls system, 2A is the energy gap and a 
is the lattice constant, which below we will put equal to 
1. The excitation energy of the SP soliton. Eg, can eas- 
ily be determined for the XK- model, i.e. if one ignores 
Sf ■ Sf^i term in the Hamiltonian (|^). Then the spin- 
Peierls Hamiltonian (J^ after the Jordan- Wigner trans- 
formation acquires a form of the Su-Schrieffer-Heeger 
Hamiltonian pO| fo r electronic Peierls materials, in which 
case Es — §A [[20|. The omitted term renormalizes the 
soliton energy, as well as the mean-field energy gap 2A, 



but these numerical changes do not play an important 
role. One should also note that the kinks are mobile 
excitations with the dispersion ^ Eg. In CuGeOa is 
estimated to be of the order of 8 lattice spacings. 

From Fig. 4 it is clear that a soliton in SP system cor- 
responds to one unpaired spin. Thus, these elementary 
excitations have ^ rather then 1 as the singlet-triplet ex- 
citations. Of course, for fixed boundary conditions the 
solitons (excited e.g. thermally or optically) always ap- 
pear in pairs. 

So far we considered the excitations in an isolated SP 
chain. Now we want to include the effects of the inter- 
chain interaction. Due to this interaction (mediated, for 
instance, by three-dimensional phonons) SP distortions 
of neighbouring chains would prefer to be phase coher- 
ent, e.g. in phase. When a kink-antikink pair of size r is 
created in one of the chains, the phase of the distortion 
between the kink and antikink is opposite to the initial 
one as well as to those on neighbouring chains, which 
would cost an energy E(r) — Zar, where a is the effec- 
tive interaction between the Peierls phases on different 
chains per one link and Z is the number of neighbour- 
ing chains (See Fig. 5a). Therefore, in the presence of the 
interchain interaction the soliton-antisoliton pair forms 
a string and Zcr may be called the string tension. The 
linear potential of the string confines the soliton motion, 
i.e. kink and antikink can not go far from each other in 
an ordered phase. 

We can use this picture to estimate the value of the 
temperature of the 3d SP transition. The concentra- 
tion of thermally excited kinks in an isolated chain is 
n = exTp{—Es/T) and the average distance between them 



is d{T) 



exp(i?s/r). At the same time, the av- 



erage distance between the kinks connected by string is 
1{T) = T/{Za). The three-dimensional phase transition 
(ordering of phases of the lattice distortions of different 
chains) occurs when l(T) d(T), i.e. 
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where we use the relation Es ^ A ^ X'J In this 
picture at T < Tsp the phases of SP distortions of dif- 
ferent chains are correlated and all solitons are paired. 
At T > Tsp local SP distortions still exist in each chain, 
but there is no long range order. Therefore, the SP tran- 
sition in this picture is of a "deconfinement" type, which 
is somewhat similar to the Kosterlitz-Thouless transition 
in 2d-systems. 

The approach described above is valid when the value 
of the interchain interaction a is much smaller than J. 
Using Eq.d) with J = lOOK and A' ^ 0.2 [| we get 



4 



a 0.04 J, which in view of the logarithmic dependence 
of Tsp on a in (^) is just enough for applicability of 
the results presented above (these are, of course, only an 
order of magnitude estimates). 



V. SOLITONS IN DOPED SYSTEMS 

As we have seen above, Zn and Si, the two most studied 
dopands of CuGeOa, lead to an effective interruption of 
spin chains into segments of finite length. The segments 
with even number of Cu ions can have a perfect SP order- 
ing, while the odd segments behave differently: one spin 
^ remains unpaired, which means that the ground state 
of an odd segment contains a soliton (similarly to what 
happens in the electronic Peierls materials |2l[]). One can 
show that the soliton is repeled by ends of the segment, 
and in an isolated odd segment the situation would look 
like in Fig. 4: the soliton carrying spin ^ would prefer 
to stay in the middle of a segment. This conclusion is 
in contrast with the usual assumption that the magnetic 
moments induced by doping are localized near the impu- 
rities. 

The situation, however, changes when we take into ac- 
count the interchain interaction. As we have seen in the 
previous section, moving a soliton along a chain costs 
an energy which grows lincary with the distance. As is 
illustrated in Fig. 5a, this provides a force pulling the soli- 
ton back to the impurity. Thus the soliton moves in a 
potential shown in Fig. 5b: it repels from the impurity 
with the potential Vim,p{r) ~ Jexp(— r/^o)i while the in- 
terchain interaction gives the potential Vconfif) ~ Zar, 
providing the restoring force. As a result, the soliton is 
located at a distance from impurity, so that, in a 

sense, we return to the traditional picture. One should 
keep in mind, however, that for a weak interchain interac- 
tion the total potential Vimp + Vconf is rather shallow and 
at finite temperature the soliton can go rather far from 
impurity. It seems that it should be possible to check 
this picture experimentally, e.g. by detailed NMR study 
of doped SP compounds (cf. the results of M. Chiba et 
al this conference). 



VI. PHASE DIAGRAM OF DOPED 
SPIN-PEIERLS SYSTEMS 

One can use this picture to describe qualitatively the 
dependence of the phase transition temperature Tsp on 
the concentration of dopands x. Similar to the treat- 
ment given in section 4, we compare an average distance 
between the kink and the nearest end of the segment 
1{T) ~ T/{Za) with the average length of the segment 
d ^ 1/x. This gives. 
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This result has two limitations. At large a;, when an 
average length of the segment becomes of the order of 
soliton size, Xjx ^ ^qj there will be no ordering even 
at T 0. Thus x ~ is an absolute limit beyond 
which the 3d ordering disappears. Using our estimate 
^0 ^ 8, such Xmax ~ 15%. On the other hand, the result 
(^ is also not valid at very small x. When an average 
size of segment d{x) ~ 1/x becomes sufficiently large, 
the thermally induced solitons become as important as 
the solitons induced by disorder. In this case the total 
concentration of solitons is 
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and one should compare 1{T) with n^J^. For a: = we 
return to Eq.(|^), while for small x we get. 



Tsp{x)^TsPmi-ax) , 



where the coefficient a is 
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One can verify these estimates more rigorously by map- 
ping the spin-Peierls system onto an effective Ising model. 
Let us associate the classical Ising variable r = ±1 with 
the two possible types of SP ordering (phases and tt), 
so that the phase (left domain in Fig. 4) corresponds to 
T = +1, while the phase tt (right domain in the same fig- 
ure) corresponds to r = — 1 . In this language a soliton is 
a domain wall in Ising variables. Since it costs an energy 
Es to create a soliton, the Hamiltonian of the intrachain 
interaction in the effective Ising model can be written as. 
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(here a is the chain index and n is the site number in 
chain). Similarly, an interchain interaction in terms of 
Ising variables has a form. 
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where the summation over 5 goes over neighbouring 
chains. One can also introduce impurities in this effective 
Ising model. The detailed treatment of this model will be 
given in a separate publication |22[ . Here we limit our- 
selves by presenting in Fig. 6 the results of the numerical 
solution of the equation for the transition temperature 
for several values of the interchain interaction a. The 
value of J was adjusted to make the transition temper- 
ature equal to 14K for each value of a. The values of a 
(from the top curve to the bottom one) are (in K) 3; 2; 1; 
0.5; 0.1. We see that the behaviour of Tsp{x) agrees with 
the (^ at large x and (^ at small x and for the values of 
the parameter a not much different from the estimates 
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made in section 4 one can obtain a reasonable form of 
the phase diagram for CuGeOa. (One should also take 
into account that each Ge is coupled to two chains, so 
that its substitution by Si introduces two interruptions 
in exchange interaction, whereas Zn interrupts only one 
chain.) 

As follows from our picture, each soliton introduced by 
doping carries an uncompensated spin i . One can easily 
show that in the vicinity of the domain wall where the 
SP order parameter is small there exist antiferromagnetic 
spin correlations (see Fig. 7). Both these correlations and 
the SP distortion change on a length scale ^o- Antifer- 
romagnetic correlations on neighbouring kinks may over- 
lap, which could, in principle, lead to the long-range an- 
tiferromagnetic ordering. Thus it is possible to obtain a 
regime in which the SP and antiferromagnetic orderings 
coexist. To study this question in detail one must also 
take into account also the interchain exchange interac- 
tion. This question is now under investigation p4|. 



VII. CONCLUDING REMARKS 

To summarize we have a rather coherent picture of 
the main properties of the SP system CuGeOs. The 
treatment given in the first part of this paper allows 
one to explain many of the features of this compound, 
which, at first glance, look rather puzzling, such as the 
strong anomalies observed in the direction perpendicu- 
lar to chains rather than parallel to them. Furthermore 
we showed how the local geometry and the side-groups 
(Ge, Si) lead to a rather detailed microscopic picture of 
the distortions in CuGeOa both above and below Tsp- 
These results are largely specific for this particular com- 
pound, although some of the conclusions {e.g. the role of 
side groups in superexchange and the importance of the 
soft bending modes) are of a more general nature. 

The results of the second part of the paper, though 
inspired by the experiments on CuGeOs, have a general 
character, e.g. the conclusions about the domain wall 
structure of the elementary excitations, confinement of 
solitons caused by the interchain interaction, disorder- 
induced solitons |23|, etc. At the same time, this general 
treatment provides a reasonable explanation of the sup- 
pression of Tgp by doping and allows to describe, at least 
qualitatively, the phase diagram of doped CuGeOa. 

We are grateful to J. Knoester, A. Lande, O. Sushkov 
and G. Sawatzky for useful comments. D. Kh. is es- 
pecially grateful to B. Biichncr for extremely useful dis- 
cussions and for informing him of many experimental re- 
sults prior to publication. This work was supported by 
the Dutch Foundation for Fundamental Studies of Matter 
(FOM). 
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FIGURE CAPTIONS 
FIG. 1. Simplified structure of Cu02 chains in CuGeOs. 

FIG. 2. Electronic orbitals of Cu, O, and Ge relevant for superexchange. 

FIG. 3. Schematic structure of Cu02 - Ge skeleton in CuGeOa. 



FIG. 4. SP soliton in the strong coupling limit (above) and in the continuum model (below). 



FIG. 5. (a) A string confining soliton to an impurity (indicated by a triangle); (b) total potential acting on the soliton 
near the impurity. 



FIG. 6. Dependence of the SP transition temperature on the concentration of dopands for several values of interchain 
coupling. 



FIG. 7. Antiferromagnetic correlations around the soliton. 
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